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Abstract. If supersymmetric particles are produced at the Large Hadron 
Collider it becomes very important not only to identify them, but also 
to determine their masses with the highest possible precision, since this 
may lead to an understanding of the SUSY-breaking mechanism and 
the physics at some higher scale. We here report on studies of how 
such mass measurements are obtained, and how the precision can be 
optimized. 
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1. INTRODUCTION 

Supersymmetry [1] has attracted a lot of attention as a natural extension 
of the Standard Model of particle physics. Not only does it solve the hierar- 
chy problem, it also has many other attractive features, the most celebrated 
of which are that it provides a natural mechanism for generating the Higgs 
potential which breaks the electroweak symmetry [2] and that it supplies a 
good candidate for cold dark matter [3]. Furthermore, if it is to be relevant 
in solving the hierarchy problem it must exhibit experimental consequences 
at the TeV-scale, and therefore can be tested by experiments at the Large 
Hadron Collider (LHC). 

If supersymmetric particles are produced at the LHC, it will become 
important to identify them and accurately measure their masses. Only an 
accurate determination of the supersymmetric particle masses and couplings 
will allow us to determine the low energy soft supersymmetry breaking pa- 
rameters. It is hoped that extrapolation of these masses and couplings to 
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high energies will provide an insight into the mechanism of supersymmetry 
breaking and, more generally, physics at the GUT scale [4]. 

Here we will discuss supersymmetric mass measurements with reference 
to one particular model of supersymmetry breaking, minimal super-gravity 
(mSUGRA) [2, 5]. In this model, the supersymmetry is broken by the in- 
teraction of new particles at high energy which are only linked to the usual 
particles by gravitational interactions; this new sector of physics is often re- 
ferred to as the hidden sector. At the GUT scale, the scalar supersymmetric 
particles are assumed to have a common mass, mo, while the gauginos have 
a common mass mi/2- The trilinear couplings are also taken to be universal 
at the GUT scale and denoted ^o- 

Furthermore, in this model R-parity is conserved. Thus, supersymmetric 
particles are produced in pairs, and the lightest one (from each decay chain) 
escapes detection. While this provides a very distinctive missing energy 
signature, it complicates the measuring of masses at the LHC since decays 
cannot be fully reconstructed. 

Instead, mass measurements rely on continuous mass distributions of de- 
cay products which attain extrema for certain configurations of the particle 
momenta that are unambiguously determined by the masses of initial, inter- 
mediate and final particles involved. These relations may often be inverted 
to give the masses of unstable particles. 



We here report on how supersymmetric mass measurements can be made 
[6, 7] by examining the mass distribution endpoints or 'edges' of the long 
decay chains (see Fig. 1) 




Figure 1: Typical cascade decay chain. 
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where g denotes the gluino, q a squark, / a slepton, whereas X2 Xi &re 
the two hghtest neutrahnos. For easier reference a subscript is given to the 
quarks and leptons when needed ("n" for "near" and "f" for "far"). 
For the decay chain (1) the following four invariants can be defined: 

mil, m,i^, m,i„ m,u, (3) 

where the subscripts arc left out when there is no ambiguity. For the longer 
chain (2) seven more invariants can be defined, but for simplicity we will 
here keep mostly to the first CcLSG ^ clS adding the gluino in any case involves 
no additional complications on the conceptual level. 

2. INVARIANT MASS DISTRIBUTIONS 
In practice, one will usually not be able to distinguish the "near" and the 
"far" leptons of Fig. 1 and Eq. (3). Therefore, two alternative distributions 
are defined: 

mqi(\ovf) = min(mg;„ , ruqi^ ) , mg,(high) = max(mqi„ , ruqi^ ) . (4) 

For the Snowmass mSUGRA benchmark point SPS la (q) [8, 6], the 
electroweak-scale masses are given as 

{rriq^, m^o, mi^, m^o} = {537.2, 176.8, 143.0, 96.1} GeV. (5) 

For these mass values, the invariant distributions for 

mil, '"^qZ(low)) '"^gZ(high)> mqll (6) 

are given in Fig. 2. For comparison, we also show the corresponding distri- 
butions for the point SPS la {(5) [6], a related mass scenario, also on the 
SPS la line [8], for which the masses are 

{niq^, m^o, mj^, m^o} = {826.3, 299.1, 221.9, 161.0} GcV, (7) 

We see that while the dilepton mass distribution has a simple triangular 
shape, and m^^^io^) for SPS la has a smooth rounded shape, the other dis- 
tributions have several abrupt changes in slope. The reason for this is that 
they are composite, with different functions describing the different inter- 
vals. In fact, these functions have recently been worked out analytically [9]. 
The great variety of these functional forms represents a complication to the 
experimental determination of the over-all endpoints. 
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Figure 2: Invariant mass distributions for the SPS la (a) [8, 6] (upper row) 
and SPS la [6] (lower row) benchmark points. 



Most of the endpoints are given by composite functions with different 
expressions for different mass constellations. For example, the endpoints of 
"^^(low) and m^^^igi^) are given by [10, 6]: 
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and the different cases given by the linear and quadratic ratios of the under- 
lying (unknown) masses: 

Case 1 : 2m? > rn?o + m-o > 2m.-am^o 
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Case 2 : m?o + rn^o > 2m? > 2m^om^o 
Xi X2 Ir ^1 ><:2 

Case 3 : m-o + "^-o > 2m^om^o > 2m~ (10) 

Xi X2 Xi X2 (fl ^ ' 

These methods for determining the masses of supersymmetric particles 
have been extensively studied in the past [11]. Our studies [6, 7, 12] ex- 
tend the earlier works by (i) discussing theoretical distributions which arise 
for different mass scenarios (examples are shown in Fig. 2), (ii) providing 
inversion formulas (by using the endpoints of the four distributions (6) the 
edge formulas can be solved for the masses m^^, m^o and m^o), (iii) 
discussing ambiguities and complications related to the composite nature of 
the endpoint expressions, (iv) extension of the method to include a gluino 
at the head of the chain. 

3. Results 

In order to estimate the precision that may be obtained at the LHC, for 
the actual masses nig, m^^, irij^, m^o and m^o, we have investigated the 
decay chains (1) and (2) for two mass scenarios, SPS la (a) of Eq. (5) and 
SPS la (/3) of Eq. (7). 

The latter has higher masses. This would be an advantage for the kine- 
matics, but higher masses also lead to lower cross sections, and thus a signal 
more easily masked by the background. The simulation was performed using 
the Monte Carlo program PYTHIA 6.2 [13], and the events were then passed 
through ATLFAST [14], a simulation of the ATLAS detector. 

A critical issue in this kind of analysis is the background. The Stan- 
dard Model background can to a large extent be removed, it turns out, 
by exploiting the fact that the two leptons must be of the same flavour. 
In most cases, the SM background events which mimic the signal, such as 
tt — > bbW^W^ hhl^l^uv produce as many different-flavour lepton pairs 
{e^jX^) as same-flavour lepton pairs (e+e", /x+jU")^. Since these two back- 
ground selections are statistically identical, the same-flavour part, which is 
collected together with the signal, can be removed statistically by subtract- 
ing the different-flavour part. Other SUSY decays will however contribute 
background events that are harder to control. The considered squark (or 
gluino) will typically be produced in association with another squark or 
gluino, which gives rise to a similar decay chain. This leads most impor- 
tantly to an additional hard jet. (We make here the usual assumption that 
the mass difference between the squarks and the sparticles into which they 

^Taus are not used because of their bad experimental resolution. 
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decay is large, producing on average the hardest jets in the event.) As it is 
not possible to know which of the two leading jets belongs to the signal, one 
will in principle select the incorrect jet in half of the cases or more. This adds 
considerably to the background and is the main reason why the experimental 
distributions tend to differ noticeably from the original distributions. 
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Figure 3: Some of the available experimental mass distributions for 
SPS la (a). The triangles mark the positions of the theoretical endpoints. 



For the two SUSY scenarios reported on here, a selection of the avail- 
able different-flavour-subtracted mass distributions ('SF-DF') are shown in 
Figs. 3 and 4. The plots are given for 300 fb~^, which is attained after three 
years at design luminosity. The Standard Model part of the samples is plot- 
ted separately in green ('SM'). Although the SM background is likely to be 
somewhat underestimated in these analyses, it is seen to have limited impact 
on the total distributions. The first four distributions of Fig. 3 correspond 
to the upper set of distributions in Fig. 2. While the experimental version of 
the niqii distribution no longer resembles its theory counterpart, appropri- 
ate cuts amend this situation for the two niqi distributions, although at the 
price of reduced statistics. The two last distributions are constructed from 
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Figure 4: Some of the available experimental mass distributions for 
SPS la (/?). The triangles mark the positions of the theoretical endpoints. 



events with two 6-tagged jets and access the gluino part of the decay chain, 
allowing for the determination of the gluino mass. 

Although there are noticeable smearing and background effects in all the 
distributions, the edges are quite pronounced and point in the close vicinity 
of the theoretical endpoint. This is also true for the SPS la {(3) distribu- 
tions, where the signal cross-section is some 25 times lower. Measuring the 
distribution endpoints is in itself a troublesome task. Apart from the statis- 
tical uncertainty, there is a (yet) largely uncontrolled systematic error from 
the signal and background hypothesis made in the actual fitting. Attempts 
at improvement are being pursued, and will be much more so if SUSY is 
discovered at the LHC. In our investigation we have therefore taken the op- 
timistic short-cut and assumed that this systematic error will be dominated 
by the controllable statistical and jet energy scale errors. 

Furthermore, rather than obtaining the SUSY masses for the specific 
endpoint measurements found from the distributions in Figs. 3 and 4 as well 
as similar distributions not shown here, these endpoints, or more precisely 
their statistical errors, were used as indicative of what the LHC will be able 
to do, and an ensemble of typical LHC experiments was generated. For each 
experiment the SUSY masses were found from a numerical fit- See [6, 10] 
for details. 

While the endpoint measurements provide us with sufficient information 
to obtain the masses, inherent to the method are two complicating effects. 
First, the endpoints relate to mass differences much more than to the masses 
themselves, hence the errors on mass differences are much smaller than on 
the masses. For the lighter masses the precision for the two quantities can 
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typically differ by around an order of magnitude. This must be considered 
a limitation of the method, since what enters further application are usually 
the masses, not mass differences. The most familiar example is the input 
to the RGE equations which on the basis of the SUSY masses allow us to 
make statements on the GUT-scale physics [4]. Another example is cross- 
section considerations. Second, there is a high probability that several sets 
of masses give the same set of endpoints. If the number of masses is equal 
to the number of measured endpoints, this is exactly true. If more endpoint 
measurements are available, the system becomes overconstraincd in favour 
of the correct solution, but due to measurement uncertainties, more than one 
solution must still often be considered. The relevant measure for comparing 
minima, is their difference in value for the mass fit. If this difference is 
small, both minima are at the same level of compatibility with the endpoint 
measurements and must be accepted. In order to choose one rather than the 
other, more information is required than what the endpoints alone provide. 
(For a discussion of these ambiguities, see [12].) 

In the case of SPS la (a) there are two competing solutions, the nominal 
one, which is in mass region (1,1)^, and a second one in region (1,2). If we 
accept solutions for which the difference in down to the global minimum 
is less than 1 (3), there is a 12% (30%) probability that the system has 
two solutions. The probability that the nominal-region solution appears is 
very high, but there is a substantial admixture of the (1,2) solution. In 
Table 1 the resulting masses for the two solution types are shown together 
with the nominal values ("Nom"). The masses are given first, then mass 
differences, with m^o selected as reference point. The better precision of the 
mass differences is apparent, especially for the lighter masses. 

In the case of SPS la (/3) three minima are competing, the nominal 
solution, a second solution with similar masses and a third solution with 
much higher masses. Again, the nominal solution is available in most cases, 
but the admixture from the other one or two is significant. The preeminence 
of mass differences over masses, in terms of precision, is once more confirmed. 
Numbers can be found in [6] . 

Both the intrinsic "weaknesses" of the endpoint method discussed above 
involve an uncertainty of the SUSY mass scale. A Linear Collider measure- 
ment of the LSP mass effectively sets this scale, so combining the measure- 
ments from the LHC with those from a Linear Collider improves drastically 
the accuracy on the SUSY masses [15]. 



^See [6] for the definition of the mass regions. 
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Table 1: SPS la (a). Nominal masses, "Nom", ensemble means, (m), and 
root-mean-square distances from the mean, a, all in GeV. 





Nom 


(1,1) 

(m) a 


(1,2) 

{m) a 


Al 

mr 

l-R 
X-2 

m~g 


96.1 
143.0 
176.8 
537.2 
595.2 


96.3 3.8 
143.2 3.8 
177.0 3.7 
537.5 6.0 
595.5 7.2 


85.3 3.4 

130.4 3.7 

165.5 3.4 
523.5 5.0 
582.5 6.8 


mr — m^o 
m^o — m^o 
mq^ - m^o 
m-g - m^o 


46.92 
80.77 
441.2 
499.1 


46.93 0.28 
80.77 0.18 
441.2 3.1 
499.2 5.6 


45.11 0.72 
80.19 0.29 
438.0 2.8 
497.0 5.4 
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